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We discuss the consequence of global polarization of the produced quarks in non-central heavy-ion 
collisions on the spin alignment of vector mesons. We show that the alignment is quite different 
for different hadronization scenarios. These results can be tested directly by measuring the vector 
mesons' alignment through angular distributions of the decay products with respect to the reaction 
plane. Such angular distributions will give rise to azimuthal anisotropy V2 of the decay products in 
the collision frame. Constraints provided by the data on the azimuthal anisotropy of hadron spectra 
at RHIC points to a quark recombination scenario of hadronization. 
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Due to the presence of a large orbital angular momentum of the parton system produced at the 
early stage of non-central heavy-ion collisions, quarks and anti-quarks are shown recently Q to be 
polarized in the direction opposite to the reaction plan which is determined by the vector of impact- 
parameter and the beam momentum. Such global quark polarization should have many observable 
, consequences such as left-right (with respect to the quarks' spin polarization) asymmetry in hadron 

spectra at large rapidity, global polarization of thermal photons, dileptons and final hadrons with 
spin. The left-right asymmetry of hadrons from the hadronization of such globally polarized quarks 
effectively gives rise to a "directed flow". However, it is difficult to disentangle it from the true 
directed flow due to interaction between produced matter and the spectator nucleons. Global hy- 
peron polarization from the hadronization of polarized quarks are predicted Q independently of 
the hadronization scenarios. Measurements of such global hyperon polarization are made possible 
by the self- analyzing power of the hyperon's parity-violating decay |2|, which produces an angular 
distribution, dN/d cos = 1 + aPjj cos 6, for its decay products with respect to the polarization 
direction. Here, a is a constant (e.g., a = 0.642 for A — > pn~) and Pjj is the hyperon's polarization. 
In practice, one needs to know not only the orientation but also the sign (or direction) of the reaction 
plane which can only be determined by the directed flow in each event class. If one sums over events 
with opposite reaction planes, the effect of the hyperon's polarization on the angular distribution of 
its decay products will cancel. While such measurements arc underway, other consequences can and 
should be studied. 

In this note, we discuss the spin alignment of vector mesons due to the global quark polarization 
in non-central A + A collisions. Similarly as in Ref. pj, we consider two colliding nuclei with the 
projectile of beam momentum p moving in the direction of the z axis. The impact parameter b is 
taken as along x, which is the transverse distance of the projectile from the target nucleus. The 
norm of the reaction plane is given by tib = p x b/\px b\ and is along y. For a non-central collision, 
the dense matter produced in the overlapped region of the collision will carry a global angular 
momentum along the direction opposite to the reaction plane (— y). Assuming that a partonic 
system is formed immediately following the initial collision, interactions among produced partons 
will lead to formation of a quark-gluon plasma (QGP) with both transverse (in x-y plane) and 
longitudinal collective motion. The global orbital angular momentum of the system will result in 
finite transverse (along x) gradient of the longitudinal flow velocity. Given the range of interaction 
Ax, two colliding partons will have relative longitudinal momentum Ap z = Axdp z /dx with orbital 
angular momentum L y ~ — AxAp z along the direction of n^. L y is estimated d] to be of the order 
of one for semi-central Au + Au collisions at yfs — 200 GeV and Aa; = 1 fm. 

Such local relative orbital angular momentum L y will lead to global quark polarization via parton 
scattering in QGP due to spin-orbital coupling. The global quark polarization via elastic scattering 
was calculated in an effective static potential model and was found 
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where E and p are the energy and momentum of the initial quark in the cm. frame of the par- 
ton scattering, p is the Debye screening mass of the quark in medium, which specifies the average 
interaction range. The polarization can be enhanced by multiple scattering. Hence, we expect a 
significant global polarization of quarks and anti-quarks before hadronization. Such global polariza- 
tion of quarks and anti-quarks will lead not only to the global hyperon polarization but also spin 
alignment of vector mesons. 

The alignment of a vector meson is described by the spin density matrix p or its element p m ,m', 
where m and m! label the spin component along the quantization axis. The diagonal elements pn, 
Poo and P-1-1 for the unit-trace matrix are the relative intensities of the meson spin component m 
to take the values 1, 0, and —1 respectively, which should be 1/3 in the unpolarized case. Since 
vector mesons usually decay strongly into two pseudo-scalar mesons, it is difficult to measure all the 
elements of p. But some of them can be determined easily by measuring the angular distributions 
of the decay products. It can be shown that, in the rest frame of V, for the decay V — ► h + h\ 
(where h and b! are two pseudo-scalar mesons), the angular distribution W(9, </>) = dN/d£l of the 
decay products is given by Q 

W[0,(j>) = -^-{coB 2 flpoo+sin 2 fl(p u +p_i_i)/2 

— sin 29 (cos 0Repio — sin 01mpio)/ V2 

+ sin20(cos0Rep_io + sin</>Imp_io)/v2 

— sin 2 #[cos(2</>)Repi_i — sin(20)Impi_i]}. 

(2) 

Here 9 is the polar angle between the direction of motion of h and the quantization axis, <f> is the 
azimuthal angle. By integrating over <j), we obtain, 

W{9) = J[(l-p O o) + (3poo-l)cos 2 0]. (3) 
Similarly, by integrating over 9, we obtain, 

W(4>) = — [l-2cos(20)Repi_i +2sin(20)Imp 1 _i]. (4) 
2n 

We see that a deviation of poo from 1/3 will lead to a non-uniform ^-distribution of the decay 
product. By measuring W(0), we can determine poo- Other elements, p w and pi_i, can be studied 
by further measuring W{6, <f>). In fact, such measurements have already been carried out in lepton 
induced reactions and hadron-hadron collisions 0, S B S H 0, Qji HH LUJ CDi LUl LL-J D3 • 

Unlike the polarization of hyperons, the spin-alignment of vector mesons, Pq , does not know 
the direction of the reaction plane since it only depends on cos 2 9 [see Eq. ©]. Therefore, one 
cannot measure the sign of the quark polarization through spin-alignment of vector mesons. On the 
other hand, one does not need to determine the direction of the reaction plane to measure the spin 
alignment which is directly related to the magnitude of the quark polarization along the orientation 
of the reaction plane. 

We now assume that quarks and anti-quarks in the QGP are polarized as described in p] and 
calculate the spin alignment of V by considering the following three different hadronization scenarios: 
(1) recombination of the polarized quarks and anti-quarks; (2) recombination of the polarized quarks 
(anti-quarks) with unpolarized anti-quarks (quarks); (3) fragmentation of polarized quarks (or anti- 
quarks). 

The picture envisaged here is the following. In a non-central A + A collision, a QGP is formed 
and the quarks and anti-quarks in the QGP are polarized. Besides them, there are also quarks 
and anti-quarks created in the accompanying processes such as the hard scattering of the partons 
and the subsequent parton cascade etc. These quarks and anti-quarks are characterized by higher 
transverse momenta and are unpolarized. Hence, there are different possibilities for hadrons to be 
produced. First, they can be produced via the recombination of the quarks and anti-quarks in 
QGP, this corresponds to the hadronization scenario (1). Second, they can also be formed via the 
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this case, we have the recombination of polarized quarks (anti-quarks) with unpolarized anti-quarks 
(quarks), and this corresponds to the hadronization scenario (2). Finally, they can also be produced 
via the fragmentation of a fast quark/anti-quark from the QGP. This corresponds to the scenario 
(3). Clearly, the three different hadronization scenarios should contribute to different kinematic 
regions. While the first scenario should play the dominant role in the low px and central rapidity 
region, the second and third should play the important roles for the intermediate px and forward 
rapidity regions respectively. 

We first consider the hadronization scenario (1) of constituent quark recombination in which both 
quarks and anti-quarks are polarized. This is likely the case for hadronization in the central rapidity 
region for low px hadrons. We take — rif, = —y as the quantization axis, and obtain the spin density 
matrix for quarks p q as, 
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and similarly for anti-quarks p q . Since the system is thermalized, there should be no intrinsic 
correlation between the quark and anti-quark in QGP. Also, since our purpose is to study the effect 
of global quark polarization, we will not go to the detail of the recombination mechanism but, just 
as people usually do |14lll5| . assume no particular correlation between the quark and the anti-quark 
that combine into a vector meson. Hence, we can calculate the spin density matrix of the vector 
meson V by making the direct product of p q and p q . After transforming it to the coupled basis, we 
obtain the normalized spin density matrix p v for vector mesons as, 
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Hence, we obtain 



V(rec) 

Poo 
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and all the non-diagonal elements are zero. Assuming P u = P c i = P fi = P c j = P q , and P s — P g , we 
obtain the results for p and K* mesons as, 



1 - P 2 

p(rec) _ q_ 

POO - 3+ P2> 



(8) 



Pm ~ 3 + P q P s ' [ ' 

We see that both p^ and p^ are smaller than 1 /3 if they are produced via recombination of similarly 
polarized quarks and anti-quarks. The non-diagonal elements are zero if there is no correlation 
between the polarization of the quark and anti-quark. 

The polarization of quark and anti-quark discussed in 1] is a low px phenomenon, since the 
polarizing interaction typically has a momentum scale of po — pLo, where 1/fi is the interaction 
range and L$ is the typical relative orbital angular momentum between two-colliding partons. When 
the initial px of a quark is much larger than po, the quark will not be polarized. But such a quark can 
still recombine with a polarized low px anti-quark to form a hadron, according to the hadronization 
scenario (2). The spin alignment for such formed vector mesons can be obtained by inserting P q = 
or Pq = into Eq. (J7J. We have then Pq Q = 1/3, even if one of the constituent quarks is polarized 
before recombination. 

Finally, we consider the hadronization scenario (3), i.e., fragmentation of a polarized quark — » 
V + X. This likely happens for quarks with large rapidities in the QGP and may play an important 
for hadrons in the forward rapidity region. The situation in this case is very much different from 
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that in scenario (1) or (2). Here, the anti-quark that combines with the initial polarized quark is 
created in the fragmentation process and may carry the information of the initial quark that induces 
this creation. This implies that the polarization of this anti-quark can be correlated to that of the 
initial quark. Since this is a non-perturbative process that cannot be calculated from pQCD, we do 
not know a priori whether such a correlation indeed exists. Fortunately, the situation here is very 
similar to e + e~ — * Z° — ► qq — > V + X, where the initial q and q are longitudinally polarized so that 
we have the fragmentation process q — > V + X. Therefore, we can compare it with the latter to 
extract some useful information. 

The 00-element of the spin density matrix for the vector mesons in e + e~ — > Z° — > qq — > V+X have 
been measured at LEP 0, 0, . The results show clearly that p^ is significantly larger than 
1/3 in the helicity frame of the vector meson (i.e. the quantization axis is taken as the polarization 
direction of the fragmenting quark) at large fractional momenta. A simple calculation 19] for Pq Q in 
e + e~ — * V + X has been carried out by building the direct product of the spin density matrix of the 
polarized leading quark (p q ) and that of the anti-quark created during the fragmentation process 
(p? as ). In the helicity frame, p q takes exactly the form as shown by Eq. J5|. The most general form 
was taken for p^ rag . The calculation is exactly the same as that for quark recombination. It also 
leads to a result of p^ for the first rank V's similar to that shown by Eq. Q. The only difference is 
that we should replace Pq in Eq. (J7J) by P^ rag , which is the polarization of the anti-quark created in 
the fragmentation process. This result has been compared with the available data0, 0, 0, 0] . It 
has been found out that, the available data can only be fitted if the anti-quark is taken as effectively 
polarized in the opposite direction as the lea ding qu ark, and the polarization is -P^ ras = — @P q , where 
[3 w 0.5 was obtained^ by fitting the dataHalllllllll- Hence, for the first rank V's, 

yjfrag) _ 1 + P P q 

For V's other than the first rank hadrons, p v = 1/3. These results can be considered as a 
parametrization of the LEP data[T3. ll5ll7lll8| . 

If the same model can be applied to the fragmentation of quarks (anti-quarks) polarized along the 
opposite direction of the reaction plane in heavy-ion collisions, then the anti-quarks (quarks) that 
are produced in the fragmentation and will combine with the leading quarks (anti-quarks) to form 
vector mesons is effectively polarized in the opposite direction as the initial quarks (anti-quarks) with 
the polarization P^ rag = —j3P q . One can then obtain a result for Pq in the same form as that shown 
by Eq. i|l(J|) . The difference is that now the quantization axis is along the opposite direction of the 
reaction plane, which is transverse to the direction of longitudinal motion. Taking the fragmentation 
of different flavors of quarks and anti-quarks into account, we obtain, for the first rank V's, 

p (frag) _ 1 + PPq , u 

pQ0 ~z~apr ( ' 



K*(f rag )_ fs l + PPj 2 n s 1 + /3P S 2 
Po ° - n s + f s 3-(3Py n s + f s 3-(3Pr [ ' Z) 

where n s and f s are the strange quark abundances relative to up or down quarks in QGP and quark 
fragmentation, respectively. Therefore, in this case of quark fragmentation, poo is always larger than 

One can measure directly the angular distribution of vector mesons' decay products with respect 
to the reaction plane and therefore determine the spin-alignment of vector mesons in non-central 
heavy- ion collisions. Such measurements will elucidate the hadronization mechanisms in the different 
kinematic regions if poo differs noticeably from 1/3. Before such data become available, one can, 
however, find constraints on poo from the measured azimuthal anisotropy of produced hadrons, 
since finite fraction of final hadrons come from vector meson decays and they have a particular 
angular distribution with respect to the reaction plane according to Eq. JSJ if poo 7^ 1/3- Such an 
angular distribution will produce an azimuthal asymmetry with respect to the reaction plane. If 
one characterizes the asvmmetrv hv the second coefficient ?>o of the Fourier transformation of the 
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angular distribution similarly as the elliptic flow study [2J}, «2 > for poo < 1/3 and «2 < if 
poo > 1/3. 

Note that the angular distribution in Eq. @ is in the rest frame of the decaying vector mesons. 
To calculate the azimuthal anisotropy of the decay products in the center of mass (cm.) frame 
of A + A collisions, one has to perform Lorentz transformation on the momentum distribution of 
the decay products from the rest frame of the vector mesons. Since quarks' polarizations should 
disappear at large px as we have argued earlier, the vector mesons' spin alignment should approach 
to poo = 1/3 at large pt- Therefore, we can assume the following ansatz for the pt dependence of 
Poo, 

Poo(pt) = P^ + (\~ Poo)- tan" 1 ^), (13) 
3 7r ao 

where Pq is the value of the spin alignment at pr = and ao sets the pr scale at which quark's 
polarization vanishes. We will use ao = 0.5 GcV/c to illustrate the effect of vector mesons' spin 
alignment on the t>2 of final hadrons. 
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FIG. 1: Azimuthal anisotropy V2 of pions from the decay of p vector mesons that have spin alignment 
according to Eq. I|13[l with po = 1/3 (solid line), (dot-dashed line) and 1 (dashed line). 



Effects of vector resonances' decay on final pions' azimuthal asymmetry has been studied recently 
[2lll22| without spin alignment. We will follow the same procedure and assume that vector mesons 
at low Pt have an exponential distribution in m-r = y/p T + m 2 with an effective temperature 
T = 200 MeV. The azimuthal anisotropy of p mesons is assumed to follow the scaling behavior of a 
recombination model 0, 

22 

1.0 + e -(p r Ao-o.35)/o.2 -0-06- (14) 

Shown in Fig.^are the azimuthal anisotropies of pions from p meson decays with three limiting cases 
of p-mesons' spin alignment. For Pq = 1/3, p-mesons are not aligned, pion distribution in the rest 
frame of the p-meson is isotropic. Therefore, V2 of pions (solid line) from the decay follows closely 
that of the p mesons. For one extreme case, Pq = 1, the angular distribution of pions prefers the 
out-plane direction and therefore V2 is negative at low px shown as dashed line. As the transverse 
momentum of the p-meson increases, the opening angle of pions from the decay in the cm. frame 
of A + A collisions becomes smaller. Eventually, V2 of pions approaches that of p-meson (solid line) 
at high pt- But it is always smaller than v-i of the p-mesons. For another extreme case, pg = 0, 
the azimuthal anisotropy of pions (dot-dashed line) is positive and larger than that of the p-mesons 
for small pt while it approaches to the p-mesons' v% at large pr from above. 

Since there are also directly produced pions, one needs to also include them in the estimate of 
the effect of p-mesons' spin alignment on the azimuthal anisotropy of the final pions. The p/7r ratio 
is measured to be 0.183 in peripheral Au + Au collisions at y/s = 200 GeV about the same value 
as in p + p at the same energy |24j . We simply assume the same value for the total number of p 

mpsnns Wf> also assnmp that ?jn for iho Hirpctlv nrnHnrerl ninns is the*, same as n-mpsnns SVinwn in 
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FIG. 2: Azimuthal anisotropy V2 of final produced pions that include both directly produced and these from 
the decay of p vector mesons that have spin alignment according to Eq. with pQ = 1/3 (solid line), 
(dot-dashed line) and 1 (dashed line). The multiplicity ratio p/V = 0.183 ,24] is assumed. The data are for 
pions from PHENIX experiment |23f . 

Fig. |21 are «2(pt) of the final produced pions with the above three limiting cases of p-mesons' spin 
alignment as compared to the experimental data [23|| which are always above the assumed v 2 of the 
p-mesons (solid line) in the scaling model 0, l22fl. The effect of vector mesons' spin alignment is 
not big due to the small value of pjix ratio. In addition to pions from p-meson decays, one also finds 
|2lL l22j ] that pions from decays of other resonances such as A and u> also enhances pions' v% at low 
Pt- For a complete analysis, one should include decays of these resonances that might also have 
similar spin alignment. However, it is clear that the data favor the scenario of Pq < 1/3 from parton 
recombination, even though the measured v 2 of pions does not provide stringent constraints on the 
value of p-mesons' spin alignment, which can be directly measured by the angular distribution of 
the decay products in the rest frame of vector mesons. 

In summary, we calculated the spin density matrix for vector mesons produced in non-central 
A + A collisions by taking into account the global polarization of quarks and anti-quarks in the 
overlapping region of the two colliding nuclei. We showed that the results for p^ depend very 
much on the hadronization mechanisms. Since mesons from the decay of aligned vector resonances 
have anisotropic angular distribution with respect to the reaction plane, they will contribute to the 
azimuthal anisotropy V2 of the final produced meson spectra. Even though experimental data on 
pions' i>2 cannot provide stringent constraints on the value of p-mesons' spin alignment, they favor 
Pq < 1/3 for the scenario of parton recombination in the central rapidity region. For the parton 
fragmentation scenario, Pq > 1/3 and vi of the decay products is negative. In this case the spin- 
alignment tends to decrease the v% of the total produced mesons. Such scenario is more likely in the 
large rapidity region and can potentially explain why the measured i>2 drops sharply with rapidity. 

It is also interesting to note that similar quark polarization effect could also happen in p + p 
collisions as pointed out in Ref. [25j and could explain the measured hyperon polarization. Therefore, 
it should also be interesting to correlate hyperon's polarization to the reaction plane in p+p collisions 
if one can similarly determine the reaction plane as in A + A collisions. However, one must also 
take into account the spin of partons inside initial beam proton as one cannot neglect them in p + p 
collisions. 
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